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Abstract 

It has by now become a standard approach to use the theory of sparse (or toric) elimination, based on 
the Newton polytope of a polynomial, in order to reveal and exploit the structure of algebraic systems. 
This talk surveys compact formulae, including older and recent results, in sparse elimination. 

We start with root bounds and juxtapose two recent formulae: a generating function of the m-Bezout 
bound and a closed-form expression for the mixed volume by means of a matrix permanent. For the 
sparse resultant, a bevy of results have established determinantal or rational formulae for a large class 
of systems, starting with Macaulay. The discriminant is closely related to the resultant but admits no 
compact formula except for very simple cases. We offer a new determinantal formula for the discriminant 
of a sparse multilinear system arising in computing Nash equilibria. We introduce an alternative notion 
of compact formula, namely the Newton polytope of the unknown polynomial. It is possible to compute it 
efficiently for sparse resultants, discriminants, as well as the implicit equation of a parameterized variety. 
This leads us to consider implicit matrix representations of geometric objects. 

Keywords. Sparse elimination; permanent; generating function; resultant matrix; discriminant; matrix 
representation 


1 Introduction 

It is today a standard approach to use the theory of sparse (or toric) elimination, based on the Newton 
polytope of a polynomial, in order to reveal and exploit the structure of algebraic systems. Consider Laurent 
polynomials in N variables: 

kf 

fi = 'y ^ C{jX i: ', C{j 0, 

3 = 1 

where x = (x\,.. . ,x/v), x e = and {o?:i, • • •,a^} C 7L N is the support of /*. Let Q * be the Newton 

polytope of fi defined as the convex hull of the support. 

A more constrained form of structure corresponds to multihomogeneous systems, which are defined on a 
product 

P ni X P” 2 X • • • X P ns (1) 
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of projective spaces: the variables are partitioned into S subsets, or blocks, so that each equation is homo¬ 
geneous in each block of n 7 + 1 homogeneous variables. In contrast, semi-mixed systems are those whose 
polynomials can be partitioned into subsets, or blocks, with fixed monomial support. In sparse elimination, 
each subset of polynomials has the same Newton polytope. In the case of semi-mixed multihomogeneous 
systems, the polynomials in each block are homogeneous of the same degree per block of variables. 

An interesting example of a semi-mixed multihomogeneous system appears in game theory and was studied 
in [20, 26]. Consider a game of S players, each with m; options or strategies. The j-th player plays a totally 
mixed strategy when the corresponding strategy is chosen randomly with nonzero probabilities p \,..., pf-f. 
A totally mixed Nash equilibrium (TMNE) is a combination of such strategies so that no player can improve 
his payoff by unilaterally choosing other (pure or mixed) strategy. We express payoff Pj of the j-th player 
choosing the i-th strategy by the polynomial 


E - 

ki,...,kj-i,kj+i,...,k s 


,U) 


(i) 




0-1) (j+1) ...J s ) 
p kj -1 Pk j+1 Pk s • 


(2) 


Here a'ff k2 ks denotes the predefined payoff of the player, assuming each player chooses the pure option 
kt G {1,..., me}. The equations imply that the payoff Pj does not depend on own strategy i, as long 
as other players do not change their strategies. Eliminating the Pj 's leads to a multilinear system with 
n,j = rrij — 1 in space (1). Actual probabilities are determined from the normalizing conditions ^”=i 'P^P = b 


for j <G {!,..., S}. 


2 Root bounds 

Bounds on the number of common roots of an algebraic system admit different types of compact formulae. 
We start with classic bounds named after Bezout, tight for (multi)homogeneous systems. The most general 
such bound is the following m-Bezout bound, stated for dehomogenized system. 

Theorem 2.1 (m-Bezout bound). Consider a system of N equations in N affine variables, partitioned into 
S subsets so that the j-th subset includes nj affine variables, and N = n\ + ■ ■ ■ + ns■ Let dij be the degree 
of the i-th equation in the j-th variable subset, for i = 1,..., N and j = 1,... ,S. If the system has finitely 
many complex roots in space (1), the coefficient of xf 1 ■ ■ ■ x^ s in 

N 

n« iiXi H-b d iS xs ) (3) 

i =1 

bounds the number of roots. For generic coefficients this bound is tight. 

We shall give a generating function for the m-Bezout number by extending results in [26, 27, 32], MacMahon’s 
Master theorem [25] is a powerful combinatorial result. 

Theorem 2.2 (MacMahon’s Master Theorem). [25] Let A = (a^) be a complex S x S matrix, x\,... ,xs 
be formal variables, and V denote the diagonal matrix with nonzero entries x\,...,xs- The coefficient of 

m 7 is - 

x l ■■ ■ x s m 

s 

JJ( a jl x l + • • • + a jS x s) nj (4) 

3 =1 
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( 5 ) 


equals the coefficient of x” 1 ■ ■ ■ x^ s in the multivariate Taylor expansion of 

f( Xl ,...,x s } = det(/ _ VA) 

around (x\ ,... , xs) = (0,0,... , 0). 


Let us consider semi-mixed multihomogeneous systems, where the number of equation blocks equals the 
number S of variable blocks, and we have a varying number of equations per block; more importantly, 
the degree per block is arbitrary. Generically, this system has a finite number of solutions. MacMahon’s 
Theorem yields a multivariate generating function for the m-Bezout bound of these systems, which could 
be helpful when one seeks root counts for a family of systems. 

Theorem 2.3. [20] Consider a multihomogeneous system on space (1) of N = n\ + • • • + ns equations, 
where the equations are partitioned into S subsets of exactly n\,...,ns equations. We assume that for 
any i,j £ {1,2,... , S'}, the polynomial equations in the i-th subset have degree a t j in the variables of the 
j-th variable subset. Let A be the S x S matrix defined by the a t j’s. Then the m-Bezout bound for the 
multihomogeneous system equals the coefficient of x" 1 • • • x^ s in the multivariate Taylor expansion of 

1/det (I-VA) 

around (xi,..., xs) = (0,0,... , 0). 


For TMNE system (2), we have an = 0, and a t j = 1 for i ^ j. Then 


I -V A = 


1 —Xl —Xl 

—x 2 1 —X 2 


-xs -X S 


-Xl 

-X 2 

-xs 1 


and 


1/det (I - V A) = 


1 - o -2 - 2cr 3 - (S - l)a s ’ 


where (for j = 2,..., S) denotes the j-th elementary symmetric polynomial in xi,... ,xg. 


Now we pass to the most general root bound in sparse elimination theory, namely mixed volume: there are 
efficient algorithms but, in certain cases, the mixed volume computation is reduced to a matrix permanent. 
The toric root bound is named after Bernstein, Khovanskii, and Kushnirenko (BKK) [1]: 

Theorem 2.4. For /i,... /tv £ Cfxf 1 ,..., x^ 1 ] with Newton polytopes Q i,..., Qn, the number of common 
isolated solutions, multiplicities counted, in the corresponding toric variety, which contains (C*) N as a dense 
subset, does not exceed MV(Qi,..., Qn), independently of the toric variety's dimension. 


It is known that mixed volume is related to the permanent. In fact, by reduction to the latter one shows 
that mixed volume is #P-complete. For a multihomogeneous system on space (1) of N = m + • • • + ns 
equations, assume the i-th equation has degree in the j-th variable block. Let A = (atj) be an IV x N 
matrix with the columns repeated nj times for each j € {1, 2,..., S}. Then the mixed volume equals [27, 
Thm 2] 


1 

n\\ ■ ■ ■ ns! 


perm A. 


A generalization to arbitrary (non-homogeneous) systems follows. For each block of variables there is a 
nj-dimensional polytope I j, each in a separate complementary space, 1 < j < S. Assume the Newton 
polytopes are direct products of scalar multiples (by a d j of the r^’s. We obtain an algebraic system on the 
product of toric varieties corresponding to Fj, of dimension rij. 
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Theorem 2.5. [20] Let A be the matrix of a- t j’s with columns repeated nj times. If there are N = ^J =1 nj 
equations, then 

S S 

MV(<2i, • • •, Qn) = perm A JJvol(Tj), Q t = JJayTj. 

j =i j=i 


3 Resultant formulae 

For the sparse resultant, a bevy of results have established matrix and, more particularly, determinantal 
formulae for a large class of systems. An introductory survey can be found in [28]. 

Resultants provide efficient ways for studying and solving polynomial systems by means of matrices whose 
determinant is a non-trivial multiple of the resultant. They are most efficiently expressed by a generically 
non-singular matrix, whose determinant is the resultant polynomial or, when this is impossible, by also 
specifying a minor M' of resultant matrix M which divides the determinant of the latter so as to yield a 
rational formula for the resultant R: 

R= det M / det M'. 

Macaulay’s classical result [24] establishes such rational formula for arbitrary dense systems, and D’Andrea, 
100 years later [7], for arbitrary systems in the context of sparse elimination. His recursive construction has 
been simplified in certain cases [17]. 

For two univariate polynomials there are matrix formulae named after Sylvester and Bezout, whose deter¬ 
minant is equal to the resultant; we refer to them as determinantal formulae. The largest family admitting 
determinantal formulae is a class of multihomogeneous systems studied in [31, 33], see also [21, Sect.13.2], 
where, each block of variables contains a single dehomogenized variable or is of linear total degree. These 
systems are moreover unmixed, in other words all polynomials have the same Newton polytope; they are 
called multigraded. 

The multigraded resultant matrices were made explicit in [10, 18], which offered determinantal formulae of 
Sylvester, Bezout and hybrid types, the latter containing blocks of both Sylvester and Bezout types. In 
[10], it is shown that there exists a determinantal pure Bezout-type resultant formula if and only if there 
exists such a Sylvester-type formula. The Bezout-type matrices generalize those identified in [6]. In [18] 
was proven the existence of determinantal formulae for systems whose Newton polytopes are scaled copies 
of one polytope. 

In [8], hybrid formulae were proposed in the mixed homogeneous case. The most complete study of hybrid 
formulae is [12] where resultants are given as a matrix determinant or a Pfaffian, i.e. a determinant square 
root, for unmixed systems of 3 dense polynomials, for up to 5 polynomials of degree up to 4, 6, or 8, and 
for 6 quadratic polynomials. 

Some studies have focused on three polynomials. In the unmixed case a determinantal formula is established 
in [23], more direct than [12] and generalizing the Macaulay-style formula obtained for the dense case in [22], 
where the numerator matrix has one row corresponding to the affine Jacobian. For Newton polygons which 
are scaled copies of a single one, the smallest Sylvester-type matrices are obtained by using a row containing 
the coefficients of the toric Jacobian [9]. In [34] they construct determinantal Sylvester-type formulae for 
unmixed systems whose Newton polygon is a rectangle from which smaller rectangles have been removed at 
the corners. 
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4 Discriminant formulae 


Discriminants are crucial in studying well-constrained algebraic systems, the system’s zero set, its singular¬ 
ities, and other properties of the system. The discriminant is closely related to the resultant but admits no 
compact formula except for very simple cases. Below we offer a new determinantal formula for a specific 
system. Recent work has shed light on the degree of mixed discriminants [5, 11] as a function of the Newton 
polytopes of the corresponding polynomials. An important open question is to obtain an explicit formula 
for multilinear systems [30]. 

In the context of sparse elimination the discriminant is defined for a set of polynomials T) with fixed supports 
At, hence it is denoted A A 1 ,...,A n (Fi, ■ ■ •, F n ) or A(Fi,..., F n ). 

Definition 4.1. The mixed discriminant A(T\,..., F n ) of n polynomials in n variables with fixed supports 
A u ..., A n in 7L n is the irreducible polynomial (with integer coprime coefficients, defined up to sign) in the 
coefficients of the Fi which vanishes whenever the system T\ = • • • = F n = 0 has a midtiple root (i.e., a 
root which is not simple) with nonzero coordinates, provided his discriminantal variety is a hypersurface; 
otherwise, the mixed discriminant is equal to the constant 1 . 

Clearly, the zero locus of the mixed discriminant is the variety of ill-posed systems, i.e., systems with at 
least one multiple root. 

We focus on systems over P 1 x P 1 x P 1 where each equation is missing one variable, motivated by the study 
of TMNE in games of 3 players with 2 pure strategies each, see polynomials (2). The bilinear system is as 
follows: 

Fi : a 0 yizi + aiyiz 0 + a 2 yoz 1 + a 4 y 0 x 0 = 0, 

F -2 : b 0 xizi + bi x\z 0 + 63 x 0 zi + b 4 x 0 z 0 = 0 , 

F 3 : c 0 X 12/1 + c 2 xiy 0 + c 3 x 0 yi + c 4 x 0 yo = 0. 

The projective coordinates (x 4 : xo), (y 1 : yo), {z\ : zq) represent the unknown probability ratios for the two 
options of each player; the coefficients a* equal the differences in the given payoffs (between the two options) 
of the first player, for the 4 combinations of pure choices of the other two players; similarly, c* equal the 
differences in the given payoffs (between their two options) of the other two players. 

The following theorem gives a 6 x 6 determinantal expression for the discriminant. The formula is of 
Sylvester-type because the matrix entries are either zero or equal to a polynomial coefficient. 

Theorem 4.2. [19] The discriminant A(T\, F 2 , F 3 ) equals 

0 0 Co ci bo b\ 

0 0 c 2 c 4 63 b 4 

co c 2 0 0 ao a 2 

det 

c 4 c 4 U U (Z 3 a 4 

bo b 3 a 0 a 3 0 0 

b\ b 4 a 2 a 4 0 0 

5 Sparse implicitization 

We argue that in sparse elimination it makes sense to consider an alternative notion of compact formula, 
namely the Newton polytope of the unknown polynomial. It typically allows computation of the polyno¬ 
mial by interpolation of the coefficients. It is possible to compute it efficiently for sparse resultants [13], 
discriminants, as well as the implicit equation of a parameterized variety. 
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This section focuses on the latter case. Given (a superset of) the Newton polytope of the implicit equa¬ 
tion, one constructs an interpolation matrix, which allows for various geometric operations to be efficiently 
computed on this matrix without need of developing the monomial representation of the implicit equation. 
The representation by interpolation matrices has been developed for parametric plane curves, surfaces, and 
hypersurfaces in the context of sparse elimination [15, 16]. It generalizes to space curves and, generally, 
objects of codimension higher than 1 by means of the Chow form, which is a generalization of resultants to 
IV-variate systems with more than N + 1 polynomials [14]. 

Another powerful matrix representation of geometric objects is based on syzygies. The corresponding theory, 
including p- bases, has been developed for parametric models [3, 4, 29]. The computation of syzygies can 
exploit the sparse elimination setting [2] and, moreover, can be achieved for point cloud models [14]. Roughly, 
syzygies of a sufficiently high degree define moving curves or surfaces that “follow” the parametric object. 
One forms the matrix expressing these moving curves and/or surfaces, which represents the geometric object 
because its rank at some point p drops if and only if p belongs to the algebraic closure of the image of the 
parameterization [3]. 
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